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5— i ■ Abstract. The recently suggested correlation between the surface temperature and 



the magnetic field in isolated neutron stars does not seem to work well for SGRs, 
AXPs and X-ray dim isolated neutron stars (XDINs; specifically the Magnificent 
Seven or M7). Instead by appealing to a Color-Flavor Locked Quark Star (CFLQS) 
we find a more natural explanation. In this picture, the heating is provided by 



magnetic flux expulsion from a crust-less superconducting quark star. Combined 
with our previous studies concerning the possibility of SGRs, AXPs, and XDINs as 
CFLQSs, this provides another piece of evidence that these objects are all related. 
Specifically, we propose that XDINs are the descendants of SGRs and AXPs. 
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1. Introduction 

Heating of neutron stars by magnetic field decay in the crust has been suggested by Pons 
et al. (2007) to explain the observed correlation between surface temperature and dipolar 
magnetic field strength of isolated neutron stars. They define a heating balance line (HBL) 
in the temperature-magnetic field diagram. Since the Soft Gamma-ray Repeaters (SGRs) 
and Anomalous X-ray Pulsars (AXPs) lie well above the HBL line and the Magnificent 
Seven (M7; a group of isolated neutron stars; see for example Haberl 2007) fall well 
below the line, here we explore an alternative explanation. Our scenario, rather than 
crustal field decay, involves magnetic flux expulsion from a superconducting crust-less 
star. The most likely compact star that can provide this is a Color-Flavor Locked Quark 
Star (CFLQS). More specifically, we employ strange quark matter in the Color-Flavor 
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Locked (CFL; Rajagopal & Wilczek 2001) phase where quarks of certain color and flavor 
pair together, resulting in a color superconducting medium. Due to the rotation of the 
star, the medium develops a vortex lattice, where the star's magnetic field is constrained 
to reside only inside these vortices (i.e. an Abrikosov lattice). 

It is generally accepted that SGRs and AXPs are the same type of objects, and it 
has been speculated before that X-ray Dim Isolated Neutron stars (XDINs) are also 
related (see Treves et al. 2000 for a review). We have previously proposed (Ouyed et al. 
2006b; Ouyed et al. 2006c; Niebergal et al. 2006) that Quark Stars in the CFL phase 
not only exist, but are manifested in the form of these three classes of astrophysical 
objects (SGRs, AXPs, and a specific group of XDINs named the M7). Using our CFLQS 
model, we present further evidence for the relation between SGRs/ AXPs and the M7 
based on an analytic prescription for the evolution of the star's effective temperature, 
T ff, and magnetic field strength. This analytic prescription is derived by considering 
vortex expulsion from the star due to spin-down from magnetic braking. We compare 
our model with observed T c g vs B and find a more natural agreement for SGRs, AXPs, 
and the M7 than the Pons et al. (2007) model. 

The paper is presented as follows: In §[2]we briefly introduce the notion of a CFLQS, 
the basis of our model, and the resulting vortex lattice formed within. We go on to 
describe the evolution of this vortex lattice in § [3] and the resulting magnetic flux ex- 
pulsion due the quantized relation between the total number of vortices and the star's 
spin-period. Lastly, in § [U we discuss XDINs (specifically the M7) and give evidence for 
their ancestral link to SGRs/ AXPs. We then conclude in §[5] 

2. CFL Quark Stars 

We assume a quark star (QS) is born with a temperature T > T c (T c is the critical temper- 
ature below which superconductivity sets in), and enters a superconducting-supcrfluid 
phase (Color-Flavor Locked; CFL phase) in the core as it cools by neutrino emission 
(Ouyed et al. 2002; Keranen et al. 2005), and contracts due to spin-down. The CFL 
front quickly expands to the entire star followed by the formation of rotationally induced 
vortices, analogous to rotating superfluid 3 He (the vortex lines are parallel to the rota- 
tion axis; Tilley&Tilley 1990). Via the Meissner effect (Meissner & Ochsenfeld 1933), the 
magnetic field is partially screened from the regions outside the vortex cores. Now the 
system will consist of, alternating regions of superconducting material with a screened 
magnetic field, and the vortices where most of the magnetic field resides. 

As discussed in Ouyed et al. (2004), this has interesting consequences on how the 
surface magnetic field adjusts to the interior field which is confined in the vortices. In 
Ouyed et al. (2006a) we performed numerical simulations of the alignment of a quark 
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star's exterior field, and, found that the physics involved was indicative of SGR/AXP 
activit\n. 



3. Magnetic Dissipation in the Crust vs Flux Expulsion 

Following the initial magnetic field alignment event is the quiescent phase, where mag- 
netic braking spins-down the QS, causing the outermost vortices to be pushed to the 
surface and expelled (Ruutu et al. 1997; Srinivasan et al. 1990). The magnetic field 
contained within these vortices is also expelled and annihilates by means of magnetic re- 
connection events near the surface of the star, causing energy release, presumably in the 
X-ray regime. The number of vortices decreases slowly as the QS spins down leading to 
continuous, quiescent, energy release that can last until the magnetic field is insufficiently 
strong to produce detectable emission. 

This scenario differs to what is expected from neutron stars, wherein the proton and 
neutron superfluids are thought to compete to push vortices to the surface, where the 
magnetic field slowly decays as it diffuses through the neutron star's crust (Konenkov 
& Geppert 2000). Pons et al. (2007) parametrize this field decay and balance it with 
blackbody cooling, thus attaining an equilibrium temperature (T e g oc B 1 / 2 ; what Pons 
et al. referred to as the heat balance line (HBL)) below which no neutron stars should be 
found. However, most of the Magnificent Seven (M7; eg. Haberl 2007) are below this pro- 
posed temperature (see Fig. [Tj), or at least require very different physical characteristics 
than their neutron star counterparts. 

In our CFLQS model, because there is only the CFL matter and no crustQ, the vortices 
are efficiently pushed to the surface where the magnetic field contained within decays 
by reconnection rather than dissipation. Thus, by balancing this heating with blackbody 
cooling we attain an equilibrium temperature proportional to B, rather than B 1 / 2 . 

This is realized by first considering spin-down due to a rotating, aligned, magnetic 
dipole, (e.g. Meszaros 1992,) 

ft B 2 R e Q 2 

ft ~ Ic? ' ( ' 

Here, ft, is the spin frequency, ft, is the spin frequency derivative with respect to time, 

B, is the magnetic field strength at the surface of the QS, R, is the radius of the QS, 

/, is the moment of inertia, and c is the speed of light. In the aligned-rotator model the 

star spins down by magnetospheric currents escaping through the light cylinder. For a 

neutron star, these currents are thought to originate in the crust. Instead, in our model, 

See simulations: www.capca.ucalgary.ca/~bniebergal/meissner/ 
2 Although it has been shown that pure CFL matter is rigorously electrically neutral 

(Rajagopal & Wilczek 2001), other work (Usov 2004 and references therein) indicates that 

a thin crust is allowed around a quark star due to surface depletion of strange quarks. In our 

model we have assumed no depletion of strange quarks, which implies a bare quark star. 
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Fig. 1. Effective temperature vs magnetic field (inferred from cyclotron resonance ob- 
servations where possible) of isolated neutron stars. The stars and diamonds represent 
SGRs and AXPs respectively, while radio-quiet, X-ray dim, isolated stars (the magnifi- 
cent seven) are shown as circles. The remaining objects are isolated radio pulsars with 
periods P > 3 s (squares) and P < 3 s (triangles). The dashed line is the heating bal- 
ance line derived by Pons et al. (2007) assuming balance between heating by magnetic 
field decay and blackbody cooling. The two solid lines are our model of vortex expulsion 
from a CFL quark star (eq. [5]). The upper and lower lines represent upper (Ai?km = 10, 
Pq = 1.5 s) and lower (Ai?km = 1, Pa — 5 s) bounds for CFL quark stars. We point out 
that the free parameters in our model are tightly constrained by observations. 

pair production from magnetic reconnection would likely supply the currents (Niebergal 
et al. 2006). 

From the quantization of angular momentum the number of vortices is proportional 
to the QS's rotation period. As given in Ouyed et al. (2004), for a sphere this relation is 
given by 

dn ~ n ' u 

where N v is the total number of vortices. 

Thus, it is easy to imagine that as the star loses rotational energy and spins down, the 
QS will lose vortices. The magnetic field contained within these vortices is also released 
from the QS. This implies that the magnetic field possessed by the QS is dependent on 
the spin period. However, the rate of spin-down is proportional to the magnetic field 
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squared (cf. Eq.[T]), so the period is also dependent on the magnetic field. Hence, the spin 
period and magnetic field are coupled, but they can be solved for independently as done 
by Niebergal et al. (2006), yielding the important relations, 

P = P (l + t/ T ) 1/3 

B = B (l + t/r)- 1/6 , (3) 

where in the above equations the subscript, 0, refers to the initial value at the time of 
the QS's birth. Also, the characteristic age, r, in units of years is calculated to be, 

where, Mqs, is the mass of the quark star, and, Pqs, its radius. 

As the magnetic field is forced outside of the star, it decays by reconnection, causing 
heating on the surface. Assuming an efficiency of 10% for the conversion to X-rays from 
reconnection, then a simple model of heating balanced with cooling gives an effective 
equilibrium temperature of the QS to be, 

kT cS ~ 1.4 x 10- 2 B 13 A< / >- 1/2 keV . (5) 

In the above equation, P13 is the magnetic field strength at the surface of the QS in units 
of 10 13 G, Ai?km is the size of the emitting region in units of kilometers, and Po is the 
initial period of the QS. 

In figure [TJ the two solid lines represent the upper (Ai?km = 10, Pq = 1.5 s) and 
lower (Ai?km = 1, Po = 5 s) bounds for the temperature given by equation ([5|). These 
bounds cover the likely extent of the parameter space, given that Ai?km is unlikely to be 
larger than the QS itself (10 km), and Po should not be much less than the SGR/AXP/M7 
period average. Thus, our QS model parameters are very physical, and tightly constrained 
by observations. The data points are from Pons et al. (2007) as is the dashed line, which 
represents their heating balance (HBL) temperature for neutron star spin-down combined 
with a two-parameter best fit for magnetic field diffusion through the crust. 

We argue that in the context of the HBL model, the M7 would require a very different 
set of parameters from other neutron stars; parameters that should be mostly uniform. 
In our QS model, the M7 share a parameter space with SGRs and AXPs, suggesting the 
two groups are the same type of objects that differ primarily in age. The other objects in 
figure [H that do not fall within the bounds of our model are regular neutron stars. Also, 
there may in fact be more objects in the gap between SGRs/ AXPs and the M7, but they 
would likely appear like regular X-ray pulsars with no persistent pulsed radio emission, 
of which there may be some unidentified candidates in the ROSAT catalogue. 
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4. The XDIN, AXP, & SGR Link 

The M7 are a class of ~ 10 6 yrs old stars possessing relatively strong magnetic field 
strengths (10 13 to 10 14 G) and exhibiting a clustering in their observed periods similar 
to that of AXPs and SGRs. Like AXPs/SGRs they show no persistent pulsed emission 
in radio wavelengths. They are also characterized by a near perfect blackbody spectrum 
(Posselt et al. 2007). The near perfect blackbody fits naturally within the framework 
of our model as the CFLQS is expected to possess no crust, but rather a bare surface. 
Moreover, the lack of radio pulsations, in our model, is a necessary consequence of the 
birth of a CFLQS, which causes the star's magnetic field to align with its rotation axis 
(Ouyed et al. 2006a). 

Although XDINs have previously been speculated to be related to AXPs and SGRs 
(see Treves et al. 2000 for a review), the most popular neutron star models were unable 
to explain period clustering and sustainment of the magnetic field. In our CFLQS model, 
after the QS's field has aligned, it will spin-down through magnetic braking as described 
in equations ([3]), and for ages on the order of ~ 10 6 yrs, we arrive at results indicative of 
the M7. As an example, if a CFLQS is born with a radius of 9 km, period of Pq — 3 s, and 
magnetic field strength of £>o = 10 14 G, then by the time it reaches ages estimated for 
XDINs it will have attained a period of 10 s and its field will have decayed to ~ 5 x 10 13 G 
(see Fig. 3 in Niebergal et al. 2006). 

Hence, by using our model with SGR/AXP parameters we arrive at M7 parameters 
after roughly less than a million years, suggesting age is the primary difference between 
SGRs/ AXPs and the M7. Also, after roughly 5 x 10 4 years the magnetic field (and 
resulting luminosity) begins to drops rapidly in our model, implying no XDIN beyond 
that age should be detectable, unless it is very close in distance (Niebergal et al. 2006). 
Recent estimates of distances to the M7 (Posselt et al. 2007) satisfy this criteria, as the 
distance varies from roughly 160 to 400 parsecs. These new distance estimates also seem 
to indicate that the number of currently observable SGRs/ AXPs is consistent with the 
seven observer] 3 ] XDINs, given their ages. 

It is worthwhile to point out that, in our model, the dipole magnetic field strength 
is given by B d = 3 x 10 19 V3PP G, which is greater than the usual field estimation for 
neutron stars by a factor of V3. Thus, vortex expulsion changes the braking index of a 
spinning dipole (n — > 4), and results in an extra factor of \/?> when predicting the star's 
magnetic field strength from its spin-period and spin-down rate. This extra factor may 
help account for the discrepancy between the M7's estimated dipole field (using the usual 
neutron star braking index; n — 3) and the observed field using cyclotron resonances (eg. 
Haberl 2007). 

3 Despite intensive searches for more objects similar to the M7, none have been found since 
2001 (Haberl 2007). 
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5. Conclusion 
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We have shown in this Letter that, in the context of CFLQS model many features of 
the Magnificent Seven can be explained and evidence for their SGR/AXP ancestry was 
presented. Specifically the evolution of: i) The effective temperature; ii) spin-period; and 
hi) the magnetic field are all predicted for SGRs/AXPs and the M7 using our CFLQS 
model. A CFLQS also has the advantage of not possessing a crust, thus its bare surface is 
able to explain the near featureless spectrum of the M7 (Pons et al. 2005). The SGR/AXP 
spectrum is naturally more convoluted as they are much further away, and so would 
likely suffer from interstellar absorption effects. Other properties of XDINs such as; (i) 
the two-component blackbody, (ii) the optical excess, and (iii) the absorption lines have 
been discussed in Ouyed et al. (2006b) in the context of CFLQS (see also Ouyed et al. 
2006c). While we have suggested possible evolutionary signatures of quark stars among 
neutron stars, signatures of their birth might have already been seen (Leahy&Ouyed 
2007). 

Acknowledgements. This research is supported by grants from the Natural Science and 
Engineering Research Council of Canada (NSERC). 

References 

Haberl, F. 2007, Astrophysics & Space Science, 308, 181 

Keranen, P., Ouyed, R., & Jaikumar, P. 2005, Astrophysical Journal, 618, 485 

Konenkov, D., & Geppert, U. 2000, Monthly Notices of the Royal Astronomical Society, 313, 66 

Leahy, D., & Ouyed, R. 2007, ArXiv e-prints, 708. larXiv:0708.1787l 

Meszaros, P. 1992, University of Chicago Press, cl992 

Meissner, W., Ochsenfeld, R., 1933, Naturwiss., 21, 787 

Niebergal, B., Ouyed, R., & Leahy, D. 2006, Astrophysical Journal Letters, 646, L17 

Ouyed, R., Elgar0y, 0, Dahle, H., & Keranen, P. 2004, A&A, 420, 1025 

Ouyed, R., Niebergal, B., Dobler, W., & Leahy, D. 2006a, Astrophysical Journal, 653, 558 

Ouyed, R., Leahy, D., & Niebergal, B. 2006b, A&A, in Press [astro-ph/0608536| 

Ouyed, R., Leahy, D., & Niebergal, B. 2006c, Submitted to A&A, |astro-ph/0611133| 

Pons, J. A., Perez-Azorm, J. F., & Miralles, J. A. 2005, Memorie della Societa Astronomica 

Italiana, 76, 518 
Pons, J. A., et al. 2007, Phys. Rev. Lett. 98, 071101 

Posselt, B., Popov, S. B., Haberl, F., Triimper, J., Turolla, R., & Neuhauser, R. 2007, 

Astrophysics & Space Science, 308, 171 
Rajagopal, K., & Wilczek, F. 2001, Phys. Rev. Lett., 86, 3492 

Ruutu, V. M., Ruohio, J. J., Krusius, M., Placais, B., Sonin, E. B., & Xu, W. 1997, PRB, 56, 
14089 

Srinivasan, G, Bhattacharya, D., Muslimov, A. G, & Tsygan, A. J. 1990, Current Science, 59, 
31 



8 Niebergal et al.: Quark Stars among compact stars 

Tilley, D. R., & Tilley, J. Superfluidity and Superconductivity, 3rd ed., Bristol, England, Hilger 
1990 

Treves, A., Turolla, R., Zane, S., & Colpi, M. 2000, The Publications of the Astronomical Society 

of the Pacific, 112, 297 
Usov, V. V. 2004, Phys. Rev. D, 70, 067301 



